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Spin-orbit (SO) splitting, iiuso, of the electron Fermi surface in two-dimensional systems man- 
ifests itself in the interaction-induced corrections to the tunneling density of states, z^(e). Namely, 
in the case of a smooth disorder, it gives rise to the satellites of a zero-bias anomaly at energies 
e — i:2ijJso- Zeeman splitting, ±ujz, in a weak parallel magnetic field causes a narrow plateau of a 
' width Se = 2u!z at the top of each sharp satellite peak. As ljz exceeds tJso, the SO satellites cross 

over to the conventional narrow maxima at e = ±2tJz with SO-induced plateaus 5e — 2a;so at the 
tops. 
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PACS numbers: 73.40. Gk, 73.23.Hk, 71.45Gm, 71.70Ej 



Introduction. A zcro-magnctic-ficld splitting |Q| of the electron spectrum in two-dimensional systems has its origin 
1— ^ in the spin-orbit (SO) coupling. The effect attracts a steady interest especially since the proposal [|| to utilize it in 
a spin-transistor device. This proposal was substantiated by the recent experimental demonstration |^,^ that the 
[ magnitude of the splitting, 2a;so: can be controlled by the gate voltage. 

' A direct consequence of a finite zero-field splitting is the beating pattern in the Shubnikov-de Haas oscillations 
D ^ These beats remain the prime tool for the experimental determination of Wso- However, the reliable results were 

^ obtained only for electrons in narrow-gap heterostructures, where the splitting is relatively strong (several meV 
^ X In the structures with small Wso (e.g., Wso ~ 0.5 meV in Ga^s-based heterostructures 0|), inferring the splitting 
d . from beating pattern becomes ambiguous, since it requires extrapolation of the beats period to a zero magnetic field. 
^ ' In addition to the beats of the Shubnikov-de Haas oscillations, a number of more delicate effects that are due 
I exclusively to the finite zero-field splitting of the conduction band have been observed experimentally Jsj-p^. Still, 
' O ■ extracting small Wgo values from the data remains a problem. 

5 ' In the present paper we identify another prominent manifestation of the zero-field splitting. Namely, we demonstrate 
that SO coupling causes anomalies in the tunneling conductance, g{V), at a, finite bias, V. The positions of anomalies 
V — ±2wso reveal directly the splitting magnitude. They emerge as satellites of the conventional zero-bias anomaly 
at V = 0. The latter is the result of electron-electron interactions modified by the disorder To the best of our 
knowledge, these SO-induced satellites represent the first case when a disorder-related effect in transport allows to 
infer the parameter of the intrinsic electronic spectrum. We note that, in the experiment, very fine features (on the 
scale of e ~ 0.1 meV p^-p^) in the tunneling density of states, i^(e), can be resolved upon the analysis of g{V) curves 
I at low temperatures T ^V. 

In two dimensions, the interaction-induced correction, Siy{e), to the free-electron density of states, vq, behaves as p^ ] 
Sh'{e) oc ln(|e|r), for e <C t^^ (diffusive regime), where r is the scattering time. This simple logarithmic form applies 
, in the presence of a tunneling electrode which causes a long-distance cut-off of the Coulomb interactions . Zeeman 
splitting of the electron spectrum in a parallel magnetic field leads to additional anomalies 6i'{e) oc ln(|e±2a;z|T) 
Remarkably, these anomalies emerge even in the weak-field limit, when ujz ^ This is because the scattering by 
the impurity potential does not mix the spin-split subbands. Thus, the wave functions of two particles with opposite 
G I spins that differ in energy by 2u!z are strictly identical. Naturally, the SO scattering from the impurities suppresses 

i the anomalies at e = ±.2llJz- 
(~| ' The situation with SO coupling is somewhat opposite. In this case, the disorder potential causes a mixing of the 
O . SO subbands. If the disorder is short-ranged, each scattering act results in the momentum transfer of the order of 
O ■ the Fermi wave vector, kp. Thus, for weak SO coupling, ujsoT *C 1, the mixing is strong, and the SO satellites of the 
^ , zero-bias anomaly are smeared out. If, however, the disorder is smooth, the momentum transfer is much smaller than 
V~j kp. In this generic regime of a small-angle scattering, the SO subbands are almost decoupled (see Fig. 1). As a 
rN result, the same reasoning that leads to the Zeeman satellites applies, so that the SO satellites are pronounced even if 
■ ujso is much smaller than . Analogously, with a smooth disorder, in the presence of a perpendicular magnetic field, 
■ ■ ■ ' electronic states separated in energy by a cyclotron quantum turn out to be strongly correlated pT[ |. Consequently, 
the interaction-induced Si'{e) is an oscillatory function of energy [^2| at low (compared to t~^) cyclotron frequencies. 
Under the condition of a small-angle scattering, the dominant contribution to dv comes from the Hartree correction 

H- 

Below we also study the evolution of the SO satellites with a parallel magnetic field. Due to a a non-trivial interplay 
between SO coupling and Zeeman splitting psl psf, this evolution is rather peculiar. Namely, weak magnetic field 
causes a fiat top of the SO satellite within a narrow interval Se — 2ujz <^ i-^so- Conversely, in a strong magnetic field, 
the Zeeman satellite p9| acquires a flat top of a width Se = 2a;so ^z- 
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Basic Equations. Microscopic origin of the SO coupling for two-dimensional electrons can be either inversion asymme- 
try of the host crystal or the confinement potential asymmetry [Ql . Unless two effects are comparable in strength 
prf , we can choose the form a{k x &)n for the SO Hamiltonian where a is the coupling constant, cr is the spin 
operator, and n is the unit vector normal to the two-dimensional plane. In a parallel magnetic field, that induces the 
Zeeman splitting 2[Jz, the Hamiltonian of a free electron has a form 

H=— V a{k X ar)n + ijj^a^, (1) 

zm 

where m is the electron mass. The spectrum of the Hamiltonian Eq. (|l|) represents two branches, so that in the 
vicinity of the Fermi surface we have 

E^{k) = EF + e^{k), (2) 

where 

e^{k) = hvF{k - kp) + fiA{k), (3) 

and 

A(fc) = ^/ + + 2wsoWz sin (f)k. (4) 

Here Vp = tikj^/m is the Fermi velocity, /x — ±1 is the branch index, losq — 2akp is the SO splitting, and 0^ is the 
azimuthal angle of fc. It is convenient to rewrite the diagonalized Hamiltonian Eq. ([l]) in the form 

i?-^ii;^(fc)A^(fe), (5) 

A. 

where the projection operators are defined as 

Kik)^U \. , MexpH^.)\ 

'^^ ' 2 V -fiexp{iipk) 1 J ^ ' 

The angle ipk is related to the azimuthal angle (j)k as follows 

tan = tan 0fc H — . (7) 

LOsO cos (pk 

In the presence of the disorder, the electron scattering time is determined by two processes, namely, intra-subband 
scattering 

1 [dp 



and inter-subband scattering 

1 f dp 



J ^Tr(\,{k)A,ip)^S{\k~p\) 5{E,{k)-EM), (8) 



j -^Tr(^A,{k)A^,{p)^S{\k-p\) S{E,{k)~E^,{p))^-L-, (9) 



where is the Fourier transform of the correlator of the random potential. Our assumption that the disorder is 

smooth can be quantitatively expressed as Xk.p ^ 1, where the parameter Xk,p is defined as 

Tr ( A,(fe)A_,(p)') = ^-^^^^''-^p) ^ {^u^^p)\ (^0) 



Correspondingly, Tr ^Ap(A;)Ap(p)^ = 1 — Xk^p is close to unity. Using Eq. (|^), the parameter Afep can be expressed 
through the angles 4>k and 4>p as 

Cj|o(tJso +C^zSin0fc)^ ^ n2 

V~ '^^^) K^k-cpp) . (11) 
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From Eqs. (^), we get the final expression for the scattering time 

11 1 m 



di^p S(kp(pp), (12) 



where we assumed that ujso, ^ -£'F- 

As it was discussed above, the satelhte anomaly in 5v{e) originates from the Hartree correction. In the case of two 
subbands, the expression for the energy-dependent part of the Hartree correction has the form 

Me) 1 a r f dq f dp f dp' P++ , wi X , 

xGf (e + u;,p + q)Gf (e + a;,p' + q)G^,{e,p)G^,{e,p'), (13) 

where V{p) is the Fourier transform of the screened electron-electron interaction. The retarded and advanced Green 
functions are the following matrices 

G«^^(.,p) = Y: ^f^l . = E KiP)G^-%,P)- (14) 

The two-particle vertex function, P , that is responsible for satellites, is determined from the standard Dyson-type 

equation with a kernel 

X (Gf(e + a;,pi +q)G^i(e,pi))S'(fcp|0p-0pj), (15) 

where the first three factors originate from the overlap integrals. With a bare spectrum Eqs. (|3|), (Q) the product 
Gf (e -I- ijJjPi -I- q)G^i{£,pi) averaged over the disorder is given by 

(Gf(e-f-tJ,pi + q)G:^i(e,pi)) = ^ 1 - i{Lu - A{p) - A{p + q))T + ihqv^T cos{(j)p - (t)g) \ (16) 

Small-angle scattering imphes that a typical q <C k^- This allows to set in Eq. ([l5| ) A(p-l-q) w A(p), Ap+q^p^+q w Ap^p^ , 
and 0p+q ~ ^p. Then the kernel Eq. (|l^) simplifies to 



K{p,p„q,u;) = f^) (1-Ap )5'(fc,|0p-0pj) _ ^^^^ 

V 27r / 1 _ _ 2A(p))T-h«7i(7UFTCOs(0p - (/)q) 

At this point we note that, upon integration over (j>p-^, the product (1 — Xp^pj^)S (^kp\(j>p — (/ipj) is proportional to 
(t^^ — T^lip)), where Tint is the inter-subband scattering time (^). Then, in the diffusive regime, (a; — 2A(p))r ^ 1, 
the solution of the Dyson equation for the vertex function Y^^{p,p' , q, uj) reads 

r++, , ^ 5(fcF|0p -(j)p'\) 

L __{p,p ,q,uj) = — ^; , (18) 

-i[uj -2A{p))t + Dq^T + T/nnt{p) 

where we have neglected a weak anisotropy of the diffusion coefficient D — v'^Ttr/2. 



In principle, the Dyson equation for P couples (weakly) this function to the other vertex functions, i.e. P 



This coupling would be important in the domain Wgo '''int- our case, uJsoTint S> 1, this coupling can be neglected. 
Substituting Eq. ( [l^ ) into Eq. (|l^) and performing integration over p and p' , we obtain 



^0 47r Jo Jo 271 Jo 2n ^i(^e - 2A{p)) + Dq^ + r^^lip) 

The fact that characteristic (j>p — 4>pi is small allows to set V{k-p\(t>p — (t>p'\) — ^(0), where V{Q) = (static 
screening). Then integration over 4>pi yields l/(mT). Finally we obtain 
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where the energy-dependent factor, £(e), is defined as 

r2iT 



Cie) 



d(f> 
2^ 



In 



(6-2A(0))^ 



(20) 



(21) 



We see that the energy cut-ofF in Eq. (21) is determined by Tj^J ^ r ^ . Note, that inter-subband scattering time Tint 
can be conveniently expressed through the conventional transport relaxation time Ttr- Using Eqs. (|^), (|ll|) we obtain 



Tint {4>) = 2Ttr 



(wso + sin ( 



(22) 



Equations (g0|)-(^2|) constitute our main result. In principle, the exchange correction to j/(e) also yields the anomaly 
at e « 2A. However, the exchange term is suppressed since it contains an extra factor A ^ t /rtr which comes from 
overlap integrals between different branches. Note also, that Eqs. (p0|)-(p2|) can be easily modified to the case when 
the crystalline anisotropy term [ p6| , (3{d'xkx — ayky), is present in the Hamiltonian (|l|) alongside with SO-term [0. 
Modification reduces to the replacement Az by Pkj,, and (p by 2(f) p^ ]. 
Analysis of the anomaly. 

(i) Zero-field limit: Uz <S ^/Ttr ^ '-^so- 
The shape of the satellite peak in Si'{e) is given by 
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(i67ri?j;^Tt^ 



In (e - 2tJso)'T^ + 



(23) 



The peak is well pronounced when LOgaTtr 3> 1. 

(ii) Intermediate fields: I /Ttr ^ ^so- 

The broadening of the SO-satellite peak is determined by the angular dependence of A((/)). The integral £(e) can be 
evaluated in this limit, yielding 



£(e) = 21n 



(|--so)r+^( 



(24) 



Remarkably, within a domain |e — 2Li;so| ^ ^zi there is a plateau in Siy(e), i.e. within this domain C{e) 

(Hi) Strong fields: l/rt,. ^ *C lOz. 
In this case, conversely, SO coupling determines the shape of the Zeeman satellite 



C{e) = 21n 



(25) 



Again a plateau at |e — 2ujz\ < Wso emerges at the top of the Zeeman satellite. 

Typical examples of the energy dependence of 6i'{(), obtained by numerical integration of Eq. (pT|), are shown in 
Fig. 2. They illustrate the successive broadening and then narrowing of the satellite peak with increasing magnetic 
field. 

Conclusion. The above consideration was restricted to the diffusive regime e ~ Wso ^ '^/t- It is known however 
p8|-p0[ , that the conventional diffusive zero-bias anomaly persists at high electron energies e ^ l/r (ballistic regime) . 
The question arises whether the SO anomaly survives in the ballistic regime. We will discuss this case qualitatively. 
Without SO coupling, the physical mechanism responsible for the formation of the zero-bias anomaly is combined 
scattering of a probe electron from an isolated impurity and a perturbation of the electron density caused by this 
impurity. The latter perturbation falls off with distance as s\Ti{2kFr) / r"^ (Friedel oscillation). Then the anomaly 
emerges as a result of the Bragg backscattering from this almost periodic potential profile. In the presence of SO 
coupling, a single impurity induces Friedel oscillation with three wave vectors |3l|| , namely 2A:f, 2k^ ± 2loso/{vw)- 
Then the wave vector for an electron with chirality +1 and energy 2a-'so above the Fermi level will satisfy the Bragg 
condition for the Friedel oscillation created by electrons with chirality —1. This is illustrated in Fig. 1. The efficiency 
of the Bragg scattering, and, hence, the anomaly at e = 2ujsoi is suppressed in the absence of magnetic field due to 
the fact that, for a given chirality, the spinors corresponding to the wave vectors k and — fc are orthogonal to each 
other. A transparent underlying physics of the ballistic zero-bias anomaly suggests that electrons with both chiralities 
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can experience Bragg scattering from the SO-specific 2fcp.-oscillation ||3l| . These resonances give rise to yet additional 
weak anomahes at energies e = iwso, that are absent in the diffusive regime. 

Note in conclusion, that experimental studies |3|] indicate that even in moderate quality Gaj4s-based samples with 
mobility ^ 2-10^ cm^/V-s (as in jl^) the typical value of the ratio r/rtr is 0.1. Then for ujso — 0.5nieV the parameter 

UJsoTtr ~ 6.6. 

As a final remark, SO coupling for two-dimensional holes is much stronger that for electrons. Therefore, the 
satellites of the zero-bias anomaly can be expected in the hole samples too. However, due to the warp of the valence 
band spectrum caused by the crystalline anisotropy, the subband splitting depends on the direction of the hole wave 
vector. This would lead to the smearing of the satellites. 
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